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A simplified view of blazars: clearing the fog around 
long-standing selection effects 
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ABSTRACT 

We propose a scenario where blazars are classified as flat-spectrum radio quasars 
(FSRQs), BL Lacs, low synchrotron, or high synchrotron peaked objects accord- 
ing to a varying mix of the Doppler boosted radiation from the jet, the emission 
from the accretion disk, the broad line region, and the light from the host galaxy. 
In this framework the peak energy of the synchrotron power (Vp ea k) in blazars is 
independent of source type and of radio luminosity. We test this new approach, 
which builds upon unified schemes, using extensive Monte Carlo simulations and 
show that it can provide simple answers to a number of long-standing issues in- 
cluding, amongst others, the different cosmological evolution of BL Lacs selected 
in the radio and X-ray bands, the larger fp eak values observed in BL Lacs, the 
fact that high synchrotron peaked blazars are always of the BL Lac type, and 
the existence of FSRQ/BL Lac transition objects. Objects so far classified as BL 
Lacs on the basis of their observed weak, or undetectable, emission lines are of 
two physically different classes: intrinsically weak lined objects, more common in 
X-ray selected samples, and heavily diluted broad lined sources, more frequent 
in radio selected samples, which explains some of the confusion in the literature. 
We also show that strong selection effects are the main cause of the diversity ob- 
served in radio and X-ray samples, and that the correlation between luminosity 
and z-'poak, that led to the proposal of the "blazar sequence", is also a selection 
effect arising from the comparison of shallow radio and X-ray surveys, and to 
the fact that high - high radio power objects have never been considered 

because their redshift is not measurable. 

Key words: BL Lacertae objects: general — quasars: emission lines — radiation 
mechanisms: non-thermal — radio continuum: galaxies — X-rays: galaxies 



1 INTRODUCTION 

Once considered rare sources, blazars, a type of radio 
loud active galactic nuclei (AGN) po inting their jets in 
the d i rection of the observer ( see e.g. iBlandford fc Reea 
119781 : lUrrv fc Padovanil 1 19951 ). are now being detected 
in increasingly larger numbers. Recent results from the 
Wilkinson Microwave Anisotropy Probe (WMAP), the 
Planck and Fermi satellites have established that blazars 
are the most common type of extra galactic sources found 
at microwave and 7-ray energies (|Giommi et all 120071 ; 
lAbdo et all l2010al : iPlanck Collaboration! 1201 J l. So far 
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abou t 3,000 blazars are known (Ma ssaro et al.1 120091 . 
1201 lh . but their number is steadily growing thanks to 
the Fermi l|Abdo et al.l l2010bl. l201lf). the opt ical Sloan 
Digital Sky Survey (SPSS: iPlotkin et alj|20icl h and the 
Planck (|Planck Collaborationll201ll ) surveys. Some faint 
blazars are also being detected as serendipi tous sources 
in Swift-XRT images (|Turrizianill20lb1 . l201lf l . 

While all blazars share the same property of emitting 
variable, non-thermal radiation across the entire electro- 
magnetic spectrum, they also display diversity. Namely, 
they come in two main subclasses, whose major differ- 
ence is in their optical properties: 1) Flat Spectrum Ra- 
dio Quasars (FSRQs), which show strong, broad emission 
lines in their optical spectrum, just like radio quiet QSOs; 
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and 2) BL Lacs, which are instead characterized by an op- 
tical spectrum, which at most shows weak emission lines, 
sometimes displays absorption features, and in some cases 
can be completely featureless. Historically, the separa- 
tion between BL Lacs and FSRQs has been made at the 
(rather arbitrary) rest-frame equivalent widt h (EW) of 
5 A (e.g. lStickel et al.lll99ll : IStocke et al.lll99ll ). However, 
no evidence for a bimodal distribution in the EW of the 
broad lines o f radio quasars has ever b een found and, on 
the contrary, IScarpa fc Falomd (119971 ") pointed out that 
radio-selected BL Lacs were, from the point of view of the 
emission line properties, very similar to FSRQs but with 
a stronger continuum. Most BL Lacs selected in the X-ray 
band, on the other han d , had very weak, if any, emission 
lines, and IStocke et al.l |l99ll ). when studying the prop- 
erties of the X-ray selected Einstein Medium Sensitivity 
Survey (EMSS) sample, had to introduce another crite- 
rion to identify BL Lacs, this time to separate them from 
galaxies. This was based on the Ca H&K break, a stel- 
lar absorption feature typically found in the spectra of 
elliptical galaxies. Given that i t s valu e in non-active el- 
lipticals is ~ 50%. IStocke et aTl (|l99lT ) chose a maximum 
value of 25% to ensure the presence of a substantial non- 
thermal continuum superposed to t he host galaxy spec- 
trum. This was later revised to 40% (|Marcha et alj|l996l ; 
lLandt. Padovani. fc Giommj|2002T ). 

Blazar classification depends then on the details of 
their appearance in the optical band where they emit a 
mix of three types of radiation: 1) a non-thermal, jet- 
related, component; 2) thermal radiation coming from 
the accretion onto the supermassive black hole and 
from the broad line region (at least in most radio- 
selected sources); 3) light from the host (giant ellipti- 
cal) galaxy. Figure \T\ represents these three components 
as red, blue and orange lines, overlaid to the spectral 
energ y distribution (SEP ) of four well-known blazars 
(from[G iommi et al.ll201ll ). The strong non-thermal ra- 
diation, the only one that spans the entire electromag- 
netic spectrum, is composed of two basic parts form- 
ing two broad humps, the low-energy one attributed to 
synchrotron radiation, and the high-energy one, usually 
thought to be due to inverse Compton radiation (see e.g. 
lAbdo et al.l [2010c). The peak of the synchrotron hump 
(i/p Cak ) can occur at different frequencies, ranging from 
about ~ 10 12 ' 5 Hz to over 10 18 Hz (see e.g. the cases of 3C 
273 or 3C 279 and MKN 501 in Fig. [J) reflecting the max- 
imum energy at wh ich particles can be accelerated (e.g. 
iGiommi et al.l [201 if ) . Blazars where Vp Cak is lower than 
10 14 Hz in their rest frame are called Low Synchrotron 
Peaked (LSP) sources, while those where 10 14 Hz < 



peak 



< 10 1 



Hz, and ^ eak 



> 10 Hz are called Inter- 



mediate and High S ynchrotron Peaked (ISP and HSP) 
sources respectively (|Abdo et al.ll2010d ). This definition 
extends the original division of B L Lacs into LBL and 
HBL sources first introduced by Pa dovani fc Giommil 

The large ^p Cak disparity between LSP and HSP 
blazars (up to five orders of magnitude) is the cause of 
very large differences between the intensity of the radia- 
tion emitted in different energy bands. For instance, for 
the same radio flux an HSP source can be a factor of a 
100 brighter in the optical band, or even a factor of a 



1,000 brighter in the X-ray band, than an LSP blazar. 
This induces very strong selection effects in blazar sam- 
ples discovered in different bands and led to some con- 
fusion when comparing, for example, the first radio and 
X-ray-selected BL Lac samples in the early 1990's. 

We feel that all the above factors, which play an im- 
portant role in blazar classification, have not been prop- 
erly taken into account so far. The purpose of this paper 
is then to propose a new hypothesis to explain the ex- 
istence of apparently different properties of the blazar 
subclasses, which solves at once many of the open issues 
of blazar research. 

Throughout this paper we use a ACDM cosmology 
with Hp = 70 km s~* Mpc" 1 , Q, m = 0.27 and Ov = 0.73 
l|Komatsu et al.ll201ll ). 



2 CURRENT STATUS: TWO TYPES OF 
BLAZAR POPULATIONS WITH WIDELY 
DIFFERENT PROPERTIES 

The two main blazar subclasses have many differences, 
which include: 

(i) different optical spectra (by definition). There are, 
however, a number of BL Lac - FSRQ transition objects, 
which include even BL Lacertae itself, the prototype of 
the class, which displays at times mod e rately strong, 
broad lines (e.g. Vcrmcule n et ail 1 19961 : ICapetti et al.l 
l2010l : iGhisellini et~al]|201ll) and 3C 279, a well-studied 
FSRQ , which can app ear nearly featureless in a bright 
state (|Pian et al.l ll999): 

(ii) different extended radio powers. Most BL Lacs 
have extended radio powers and morphologies consis- 
tent with those of Fanaroff-Riley (F R) type I, while 
basic ally all FSRQs are FR Il-like l|Urrv fc Padovanil 
1995], and references therein). However, despite the dif- 
ficulty of classifying the radio morphology of sources 
with their jets forming a small angle with respect to the 
line of sight, some radio-selected BL Lacs are known to 
posses an FR H-like st ructure (e.g. fKollgaar d et al.ll 19921 : 
iRector fc Stockell200ll ): 

(iii) very different redshift distributions. FSRQs, simi- 
larly to radio quiet QSOs, are typically found at redshifts 
~ 1 — 2, and up to ~ 5.5, while BL Lacs are usually much 
closer with very few cases at z ^ 0.6. This is graphically 
shown in Fig. [2] which plots the redshift distribution of 
all the FS RQs and BL Lacs inc l uded in the third edition 
of BZCAT (|Massaro et alj|2009l.l201ll ). the largest compi- 
lation of blazars currently available. It must be stressed, 
however, that a large fraction of BL Lacs (~ 43% in BZ- 
CAT and > 50 - 60% of the BL Lacs in the Fer mi 1 and 
2 year AGN catalogs: lAbdo et al.|[2010al . 12011b?) have no 
measured redshift, due to the lack of any detectable fea- 
ture in their optical spectrum, despite the use of 10-m 
class optic al telescopes for th e spect roscopy identification 
campaign l|Shaw et al.ll 20091 . l201Ch ; 

(iv) different cosmological evolutions. Detailed analy- 
ses of both radio selected and X-ray selected samples 
have led to the conclusion that the cosmological evo- 
lution of the two blazar subclasses are very different, 
with FSRQs evolving strongly (again similarly to ra- 
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Figure 1. The SEDs of four representative blazars: two FSRQs, 3C 273 and 3C 279, and two BL Lacs, MKN 501 and BL Lac. 
The lines in color denote the three main components of blazars SEDs, namely non-thermal radiation f rom the jet (red), emission 
from the disk and from the broad line region represented by the composite QSO opti cal spectrum of Vanden Berk et alj ll200lh 
(blue), and light from the host galaxy, represented by the giant elliptical template of iMannucci et ah I ll200ll ) (orange). The two 
vertical lines indicate the optical observing window (3800 — 8000 A). See text for details. 




Figure 2. The redshift distribution of FSRQs (1676 objects, 
solid line) and BL Lacs (537 objects , dashed line) in the thir d 
edition of the BZCAT catalogue jMassaro et alJ|200Sl [201 ll) . 
The two distributions are obviously different at the > 99.99 
level, with means of 1.4 and 0.36 respectively. About 400 BL 
Lacs with no redshift determination in the BZCAT catalogue 
have been excluded. 



dio quiet QSOs) and BL Lacs evolving at a similar, or 
perhaps lower rate, in the radio band, or even show- 
ing no or n e gative evolution in the X-ray band (e.g . 



Stickel et al 



Rector et al 



1991]; iGiommi, Menna fc Padovanl Il999l; 
2000;' lBeckmann et al.ll2003l ; |Padovani et all 



20071 : lG iommi et al.ll2009l . and references therein); 



(v) widely different mix of FSRQs and BL Lacs in ra- 
dio and X-ray selected samples, with the latter typically 
including a much larger fraction of BL Lacs than the for- 
mer. In fact, while only ~ 15% of WMAP5 blazars are 
BL Lacs (Section I4.1|l . this fraction is instead ~ 70% in 
the EMSS com plete sample, which includes 41 BL Lac s 
and 15 FSRQs ^Rector et al.ll200d : IPadovani et alfeOOof ) : 

(vi) widely different distributions of the synchrotron 



energy Vpc^- The rest-frame 
of FSRQs is strongly peaked at 

,.S \ _ 1r ,13.1i 

peak 

ues 



peak 

low 



distribution 
energies ( ( 



10 1 



Hz) and never reaches very high val- 
peak ~ 10 14 ' 5 Hz) independently of the selection 
method l|Giommi et al. l201ll ). while the ^p ea k distribu- 
tion of BL Lacs is shifted to higher values by at least 
one order of magnitude. It can also reach values as high 



peak 



10 Hz and its shape varies strongly depend- 
ing on the selection band (that is radi o, X-ray or 7-ray: 
lAbdo et al.ll2010d : IGiommi et alll201ll ) . 
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Some of these differences have been explained by so- 
called unified schemes, which posit that BL Lacs and FS- 
RQs are simply FR I and FR 11 radio galaxies with their 
j ets forming a small ang le with respect to the line of sight 
(|Urrv fc Padovanlll995r ). (Radio galaxies would then be 
the "parent" population of blazars). Due to relativis- 
tic beaming this has enormous effects on their apparent 
emitted power and luminosity functions (LFs) and can 
explain their different extended radio powers and, partly, 
their cosmological evolutions. However, unified schemes 
per se cannot account for transition objects, and the dif- 
ferent evolution of radio and X-ray selected BL Lacs and 
"peak distributions. Our new hypothesis builds upon uni- 
fied schemes by adding dilution and strong selection ef- 
fects as new, vital components. 



3 THE EW OF FSRQS AND RADIO QUIET 
QSOS 

The existence of BL Lac - FSRQ transition objects and 
the recent finding that unresolved radio quasars (i.e. 
FSRQs) appear to be redder than radio quiet AGN 
l|Kimball et al.l201ll ) suggest that the strong non-thermal 
jet emission in blazars may have a significant impact on 
the shape of their blue bump and on the EW of their 
broad lines. In particular, given this extra continuum 
component, the EWs of broad-lined blazars (FSRQs) 
should be systematically lower than those of radio-quiet 
AGN. 

To test this, we extracted EW data for samples 
of radio-quiet QSOs and blazars from the SDSS DR- 
7 spectral database. In order to have a sample that is 
not excessively large but still representative of the en- 
tire data set, we considered only radio-quiet QSOs with 
28° 6n ^ 30°, thus limiting the sample size to about 
650 sources. To obtain the EW of a sizable sample of 
blazars we considered all FSRQs in BZCAT with 1.4 GHz 
flux > 300 mjy (to simulate a radio flux-limited sample) 
and in the WMAP5 catalogue. For each object in the 
three samples we retrieved, from the SDSS on-line sys- 
tem, the line measurements parameters (e.g. EW, xt °f 
the fit etc.) of the strongest emission lines in quasar spec- 
tra, that is Lya, C IV, C III, Mg II, H/3, Ha. Details on 
the fitting algorithm can be found at the SDSS websitfl 
To avoid problematic cases, for our analysis we consid- 
ered only SDSS fits to emission lines with xt ^ 2.5. 

Figure [3] shows the distribution of the EW of the 
Mgll emission line, for which we have better statistics, 
for the sample of radio quiet SDSS QSOs (663 objects), 
BZCAT (126 objects) and WMAP5 (35 objects) FSRQs. 
The EW of both blazar samples is smaller than that of 
radio quiet QSOs, with median values respectively equal 
to 16.4 (BZCAT), 14 (WMAP5), and 18.4 A (radio-quiet 
QSOs). A Kolmogorov-Smirnov (KS) test gives a prob- 
ability < 0.1% that both FSRQ samples have the same 
distribution as that of radio-quiet AGN. Moreover, we 
found that this dilution increases further (i.e. the median 
EW gets smaller) if we make a cut at higher radio fluxes 
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Figure 3. The distribution of the Mg II rest-frame EW for 
a sample of 663 radio quiet quasars (solid histogram) com- 
pared to the distributions of the Mg II EW of 35 WMAP5 
and 126 FSRQs in BZCAT (dashed and dotted histograms 
respectively) with SDSS-DR7 spectrum and xl ^ 2.5. The 
blazar EW distribution is significantly shifted towards lower 
values compared to radio quiet QSOs. See text for details. 



in the BZCAT FSRQs. Similar results were obtained 
for the other emission lines. This demonstrates that the 
broad lines of blazars are significantly diluted compared 
to those of radio quiet QSOs. The mean absolute magni- 
tudes of the blazar and radio-quiet samples agree within 
~ 0.2 magnitudes, which means that the difference can- 
not be due to the well-known anti-correlation between 
EW a nd absolute m agnitude, the so-called "Baldwin ef- 
fect" (|Baldwinlll977h . 

We note that a number of papers have stated that 
the optical spectra of radio-quie t and radio-loud Q SOs 
are not significantly different (e.g. lFrancis et al Il993l . and 
references therein). However, the radio-loud samples on 
which this conclusion was based were still optically se- 
lected and required the presence of a UV excess and/or 
strong emission lines. Therefore, these radio-loud samples 
were biased towards objects with strong blue bumps. 



4 A NEW, MUCH SIMPLER, SCENARIO 
FOR BLAZARS 

The considerations made in the previous sections lead 
us to propose a new, very simple scenario where the ob- 
served blazar optical spectrum is the result of a combi- 
nation of an intrinsic EW distribution and the effects of 
three components: a non-thermal, jet related one, a ther- 
mal one due to the accretion disk, and emission from the 
host galaxy. Different mixes of these components deter- 
mine the appearance of the optical spectrum and there- 
fore the classification of sources in FSRQs (dominated 
by strong lines), BL Lacs (with diluted, weak lines, if 
a standard accretion disk is present), and radio-galaxies 
(where the host galaxy swamps both the thermal and 
non-thermal nuclear emission present in blazars). The 
other novel component is a single LF whose evolution 
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depends on radio power. This idea was tested through 
extensive Monte Carlo simulations, as described below. 



4.1 Simulation ingredients 

Our simulations include the following ingredients, which 
we kept as simple as possible and tied as much as possible 
to observational data: 

(i) Luminosity function We derive the LF and 
evolution of blazars at 41 GHz from the Wilkin- 
son Microwave An isotropy Probe (WMAP5) sample 
(| Wright et alj|2009h . Given that the mean blazar radio 
- sub-mm spectral index is ~ jGiommi et al.ll2009I ) this 
is practically equivale nt to a radio-se l ected sample. We 
extend on the work of iGiommi et all (|2009l ), which was 
based on WMAP3, and define a flux-limited sample of 
high Galactic latitude sources (fnoHz zZ 0.9 Jy, |6u| > 
15°) including 161 FSRQs, 29 BL Lacs, and 10 blazars of 
unknown type. By applying a maximum likelihood tech- 
nique to the WMAP5 blazars (see, e.g. iPadovani et al.l 
l201ll . for details) we obtain, together with the evolution 
discussed below, a best-fit local LF <f>(P) cx P -3 (in units 
of Gpc~ 3 P- 1 ) between 1.9 x 10 24 and 4.2 x 10 27 W/Hz, 
which we assume in our radio simulations. 

(ii) Cosmological evolution Powerful (P r > 10 26 
W/Hz) radio sources are known to display a strong, pos- 
itive evolution at moderately l ow redshifts follo wed by a 
decline at higher redshifts (e.g. IWall et al1l2005l . and ref- 
erences therein). We parametrize this behaviour with a 
simple model of the type P(z) = (l+z) k+l3z , which allows 
for a maximum in the luminosity ev olution followed by a 
declin e. This was first suggested b y Wall, Pope fc Scottl 
(|2008h and applied bv lAiello etahl (|2009h to a sample of 
SwiftfBAT blazars. A maximum likelihood technique ap- 
plied to the WMAP5 blazar sample allows us to derive 
k — 7.3 and ft = —1.5 in the — 3.4 redshift range (which 
implies a peak at z ~ 1.85), which we assume in our sim- 
ulations. The case of pure luminosity evolution (/? = 0) 
is excluded with very high significance (P > 99.99%). 



Lower luminosity (P r < 10 W/Hz), mostly FR I radio 
sources are known to display a much weaker cosmologi- 
cal e volution, which reaches ~ zero at P r < 10 25 W/Hz 
(e.g. iGendre. Best fc Wall 2010], and references therein). 
We took this into account by using the radio LFs of BL 
Lacs and FSRQs derived from those of FR Is and FR 
lis an d based on the beaming model of lUrrv fc Padovanil 
(|1995|K which agree well with those of recent b lazar sam- 
ples (jPadovani et al.l[2007l : lGiommi et al.ll2009r i. We then 
used the fraction of beamed FR I blazars in bins of ra- 
dio power to simulate the fraction of non-evolving radio 
sources as a function of power. This fraction is equal to 
1 for P r ^ 5 x 10 24 W/Hz, decreases monotonically with 
power, and reaches for P r ^ 5 x 10 27 W/Hz. 

(iii) Non-thermal component To represent the 
non-thermal/jet component, we assume a simple homo- 
geneous synchrotron s elf-Compton model (SSC, see, e.g. 
iTramacere et al1l2009l . and references therein) with rela- 
tivistic electrons distributed as a power l aw at low ener- 
gies and as a log-parabola at high energies (|Massaro et al.l 
I2004L [20061. This model represents well the synchrotron 
part of the observed SEDs, which always extends at least 




Figure 4. The distribution of the Lorentz factors of the elec- 
trons radiating at the peak of the synchrotron SED used for 
the simulation, which also assumes a magnetic field of B=0.15 
Gauss and a gaussian distribution of Doppler factors with 
(<5)=15. 



to the optical band where the classification of blazars as 
FSRQs or BL Lacs occurs. As for the inverse Compton 
emission, which can be important in the soft X-ray band, 
we set the Compton dominance so as to reproduce the 
observed / x //r in FSRQs. This is sufficient for our pur- 
poses, since, complicating the emission with additional 
components, like thermal emission from accretion or in- 
verse Compton on an external field of photons, would 
only modify somewhat the amount of observed soft X- 
rays in LSP blazars but would not change the composi- 
tion of the samples, nor alter any of our conclusions. 

The Lorentz factors of the electrons radiating at the 
peak of the synchrotron SED component (7 pea k) are 
distributed as shown in Fig. [4] The range of 7 poa k (~ 
2.5 — 4.5) is that expected for typi cal paramete r s of th e 
SSC model as shown in Fig. 36 of lAbdo et ail (|2010d ). 
The particular shape of the distribution was chosen so as 
to empirically reproduce the observed v^ eA distributions 
in radio and X-ray selected samples of blazars. The skew- 
ness to lower values is similar to that of the ,f x //r (a proxy 
for t^p pai,) distribution adopted by IPadovani fc Giommil 
(|1995r i to unify X-ray selected and radio-selected BL 
Lacs. As regards the Doppler factor, we assumed a mean 
value of 15, which was chosen to be consistent with 
the mean s u perlu minal speed /3 app ~ 12 obtained by 
iLister et al.l (|2009l) , and with the typi cal Lorentz factor 
r ~ 15 derived bv lHovatta et all (|2009T ) (since for the an- 
gle that maximizes the apparent velocity S ~ /3 apP ~ T). 
(iv) Accretion Disk and broad emission lines 
We use the quas ar spectral template of 
IVanden Berk et al.l (|200ll ). which has been built us- 
ing an homo geneous dataset of ov er 2,200 SDSS spectra 
(see Fig. 4 of lGiommi et aUbOTll ). A standard accretion 
disk is likely to be present only in so-called "high 
excitation" radio galaxies (HERGs), while it appears 
not to be there, or be less efficient, in low-excitation 
ones (LERGs). Almost all FR Is are LERGs, while 
most FR lis are HERGs, although there is a population 
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of FR II LE R Gs a s well. IChiaberge et al.1 (|l999l ) and 
iDonato et alJ |2004) have in fact suggested that the 
obscuring torus, which is required by AGN unification 
schemes, is absent in FR Is, based on the high optical 
core detection rate and low X-ray intrinsic absorption, 
respectively. This would mean no accretion disk as well, 
otherwise one would see broad lines in their optical 
spectra, while, apart from a handful of objects (e.g. 
3C 120), ba s ically no FR I displays broad optical lines. 
I Evans et al.l l|2006l ) have shown that the accretion flow 
luminosities of FR lis are typically several orders of 
magnitude higher than those of FR Is. It then looks like 
LERGs, which means all FR Is and some FR lis, either 
do not possess an accretion disk, or if the disk is present 
is much less efficient than in FR lis (i.e. of the Advection 
Dominated Accretion Flow [ADAF] type). We have then 
associated the presence of a standard accretion disk only 
with beamed FR II sources and assumed that all those 
with an FR I parent (the non-evolving sources) have no 
disk (but see Section [6] for an alternative scenario). 

(v) Equivalent width distributions The intrinsic 
(before dilution) distributions of the EW of the broad 
lines (Lya, C IV, C III, Mg II, H/3, Ha) have been as- 
sumed to be those of the radio quiet QSOs included in the 
SDSS database described in Section[3] We assumed Gaus- 
sian distributions characterized by the measured means 
({EWhcH 200 A, (EW H/ 3>= 23 A, (EW Mg -n)= 18 A, 
(EWc-ra)= 16 A, (EWc-iv)= 20 A, (EW Lya )= 47 A) 
and dispersions for the various lines. 

(vi) The disk to jet power ratio The disk and 
jet components in blazars are known to be corre- 
lated and possibly of the same order of magnitude, 
although there are u ncertainties associated with esti- 
mating them (see, e. g. iD'Elia. Padovani. fc Landtll2003l : 
iGhisellini et~al]l201ll ). We are interested in the somewhat 
simpler question of determining how the luminosity of 
the accretion disk (blue bump intensity at 5000 A) scales 
with radio power (at 5 GHz). The relevant data were de- 
rived by using the very large amount of multi-frequency 
information included in public databases and the tools 
that are now available to analyze SEDs (|Stratta et al.l 
l201ll ). Figure [1] gives some examples of representative 
objects. The amount of thermal flux in each FSRQ was 
estimate d by matching the compos ite optical QSO spec- 
trum of IVanden Berk et all (|200ll ) to the SED data in 
the blazar rest-frame. We have done that through a care- 
ful visual inspection of each SED and by adjusting the 
intensity of the composite QSO spectrum until it over- 
lapped well to the observed emission. The use of this 
manual approach was necessary as the heterogeneity of 
the available data and flux variability does not allow the 
implementation of a robust automatic procedure. Fig. [1] 
gives examples of the matching of the composite QSO 
spectrum (blue line) to the data for the case of 3C273 or 
3C279. In those objects where the available optical/UV 
data are limited to a magnitude in one or two colors, we 
matched the composite QSO spectrum to the flux level 
corresponding to the available magnitude(s), taking into 
account of redshift. 

An upper limit was instead estimated for BL Lac ob- 
jects by placing the composite QSO spectrum in the SED 
at an intensity such that the optical lines would not 
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Figure 5. ct r -BlueBump> the slope between the radio lumi- 
nosity at 5 GHz (non-thermal boosted component) and the 
optical luminosity at 5000 A due to the blue bump (thermal 
non-boosted component estimated from the blazars' SEDs as 
described in the text), is plotted against radio luminosity for a 
large number of blazars from several radio and X-ray selected 
samples. The solid line is the linear relationship used for our 
simulations (see text for details). 



be detectable in the optical spectrum (typically a fac- 
tor ten below the observed flux. Fig. [1] illustrates the 
case of MKN501). This was done for a large number of 
blazars selected in four surveys, two radio-selected (Deep 
X-ray Radio Blazar Sample [DXRBS] and WMAP 5) 
and two X-ray selected (EMSS and Swift-B AT). The 
results are plotted in Fig. [5] where we can see that 
otr-BiueBump (the slope between the 5 GHz luminosity 
and the blue bump intensity at 5000 A and defined by 
£disk = £r(^ 5000 ^/^5Gffz)~ ar ~ BluoBump ) correlates with 
radio luminosity, although with a large scatter. We in- 
terpret this as the result of relativistic beaming, in the 
sense that the larger the radio luminosity (and therefore 
the beaming amplification) the larger the ratio between 
non-thermal and thermal radiation. 

We adopted a simple linear relationship between the 
two variables (a r -BiueBump = 0.04*log(L ra dio) — 0.39, see 
Fig. O and assumed a Gaussian distribution around it 
with a dispersion of 0.1. This has been derived by fit- 
ting the WMAP5 and EMSS points in the Figure (see 
below). We did not take into account the lower limits on 
(Xr-BiueBump since the large majority of them are at low 
radio power ( 10 26 W/Hz) and in our scenario these 
sources do not have a standard accretion disk. 

(vii) Host gal axy 

Following IScarpa et all |2000l ) and lUrrv et al.l l|2000t ) 
we have assumed that the host galaxy of blazars is a giant 
elliptical with fi xed absolute magnitude of Mr = — 22.9 
as estimated bv lSbarufatti. Treves, fc Falomol l|2005l ) us- 
ing the same ACDM cosmology adopted in this paper. 
These authors have shown that the dispersion around 
this value is less than 1 magnitude and therefore such 
a giant elliptical galaxy can be considered a standard 
candle. For the spectral shape we used the galaxy tem- 
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plate of lMannucci et al. ( 2001), who derived it combining 
the data from 28 local elliptical galaxies observed in the 
wavelength range 0.12 — 2.4 /jm. 

Figure [TJ shows this template superposed to the SED 
of four well-known blazars. 

4.2 Simulation steps 

Our Monte Carlo simulations consist of the sequential 
execution of the following steps: 

(i) draw a value for the radio luminosity and redshift 
of a simulated blazar based on the luminosity function 
and evolution described above; 

(ii) draw a value of the Lorentz factor of the electron 
radiating at the peak of the synchrotron power (7 pea k) 
from the distribution shown in Fig. [4] which is assumed 
to be independent of blazar luminosity; 

(iii) Calculate the peak of the synchrotron power in the 
source rest frame. In our simulations we assume a simple 
SSC model, and therefore v V£a k = 3.2 x 10 6 7p eak B<5. The 
value of 7 P eak is derived in step ii), the magnetic field 
is fixed to B = 0.15 Gauss, and the Doppler factor 8 is 
randomly drawn from a gaussian distribution with (S) — 
15 and a= 2. 

(iv) calculate the observed radio flux density (from the 
radio luminosity and redshift from step i)) and the non- 
thermal emission in the optical and X-ray bands under 
the assumption that the spectral shape of the observed 
emission is a log parabola around Vp ell k and that the low 
energy part of the SED (cm and mm wavelengths) is a 
power law, as is typically seen in bl azars (see e.g. F ig.[T]for 
some representative examples, or iGiommi et all (|201lh 
for a much larger sample of blazar SEDs) ; 

(v) accept the source if its flux in the band under con- 
sideration (radio or X-ray) is above the flux limit chosen 
for the simulated survey; 

(vi) add an accretion (blue bump) component as de- 
scribed above (only for beamed FR II sources), re-scaling 
the SDSS quasar template to this value; 

(vii) draw a value of the equivalent width of Lya, C IV, 
C III, Mg II, H/3, Ha starting from the EW distribution 
observed in the SDSS radio quiet QSOs; 

(viii) add the optical light of the host galaxy assuming 
a standard giant elliptical observed in blazars; 

(ix) calculate the total optical light and the observed 
equivalent width of all the broad lines considered by tak- 
ing into account the dilution due to the non-thermal and 
host galaxy optical light; 

(x) classify the source as an FSRQ if the rest-frame 
EW of at least one of the broad lines that enter the 
optical band in the observer frame (which we assume 
to cover the 3, 800 - 8, 000 A range) is > 5 A. Oth- 
erwise, the object is classified as a BL Lac, unless the 
host galaxy dominates the optical l ight causing the Ca 
HfcK break to be larger than 0.4 (|Marcha et all Il996l ; 
lLandt. Padovani. fc Giommil I2002T ). in which case the 
source is classified as a radio galaxy. A BL Lac whose 
maximum EW is < 2 A, or the non-thermal light is 
at least a factor 10 large r then that of the host galaxy 
(|Piranomonte et al.| [2007). is deemed to have a redshift 
which cannot be typically measured. 



It is important to stress that the scope of our simu- 
lations is not to reproduce all the observational details. 
While that could be possible in theory, in practice it 
would require a large number of parameters and some 
speculations. Our approach is instead to keep the number 
of assumptions to a minimum, with the aim of allowing 
us to obtain robust, almost model-independent conclu- 
sions. Our main results are in fact independent of the 
simulation details (see Section \5A\i . 

4.3 Simulations of radio and X-ray surveys 

We simulated a radio flux density limited survey with 
/ ^ 0.9 Jy, to match the WMAP5 sample, and an X-ray 
flux limited survey down to 5 x 10~ 13 erg cm' 2 s _1 in 
the 0.3 — 3.5 keV band, in order to be able to compare 
it with the EMSS. To ensure good statistics each simula- 
tion run included 10,000 sources. In the X-ray case, since 
radio powers reach lower values than in the radio case, 
we extrapolated the radio LF down to 1.9 x 10 23 W/Hz 
assuming the same slope. 



5 COMPARING SIMULATIONS AND REAL 
DATA 

In this section we make a detailed comparison between 
the results of our simulations and the observational data 
in terms of fractions of BL Lacs and FSRQs, redshift and 
ff„i, distributions, and cosmological evolution using the 
V/Vm test (I Schmidt! 1 1 968T ) . where V is the volume out 
to the source and V m is the the volume at the distance 
where the object would be at the flux limit. We also show 
where our simulate d sources end up on the plot used by 
iFossati et al.l (|l998l ) when they first described the "blazar 
sequence" . 

5.1 Radio flux density limited survey 

Table [T] summarizes our main results by giving the num- 
ber of sources per class, their mean redshift, and (V/Vm)- 
The number in parenthesis refers to the BL Lacs with 
measurable redshift, to which the mean redshift and 
(V/Vm) pertain. About 3/4 of our sources are classified 
as FSRQs, with the fraction of BL Lacs being ~ 19.8% of 
blazars, which is consistent with the value of 15.3l3Q% 
in the WMAP5 sam ple (where th e la errors are based on 
binomial statistics: lGehreldfl986h . A small fraction (5%) 
of the simulated blazars are classified as radio galaxies. 
These are bona-fide blazars misclassified because their 
non-thermal radiation is not strong enough to dilute the 
host galaxy component. 

The mean redshift for our simulated FSRQs is in 
good agreement with the WMAP5 value of 1.13, while 
for BL Lacs this is slightly larger than the WMAP5 mean 
(0.55). 

FigE] shows the overall good agreement between our 
simulated redshift distributions (where we have only in- 
cluded sources with a measurable redshift) and the ob- 
served ones. 

63% of our BL Lacs have a redshift determination, in 
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Figure 6. Top panel: the redshift distribution of the WMAP5 
FSRQs (solid histogram) compared to that of FSRQs in a sim- 
ulation of a radio survey (dashed histogram). Bottom panel: 
the redshift distribution of the WMAP5 BL Lacs (solid his- 
togram) compared to that of BL Lacs in a simulation of a 
radio survey (dashed histogram) 



Table 1. Results from a simulation of a radio flux density 
limited survey (0.9 Jy) 



Source type 



Number of 





sources 


(*) 


{V/V m ) 


FSRQs 


7,587 


1.24 


0.64 


BL Lacs 


1,879 (1,191) 


0.87 


0.60 


Radio galaxies 


534 


0.04 


0.48 


Total 


10,000 


1.13 


0.63 



excellent agreement with the WMAP5 value of 69_ 20 %- 
79% of the BL Lacs (68% of those with redshift) have a 
standard accretion disk and are therefore broad-lined but 
are classified as BL Lacs only because their observable 
emission lines are swamped by the non-thermal contin- 
uum. 

As regards (V/V m ), our simulated mean values agree 
with the observed ones of 0.62 ± 0.02 and 0.63 ± 0.05 
for WMAP5 FSRQs and BL Lacs respectively (cf. also 
the value of 0.60 ± 0.05 for the 1 Jy BL Lac sample: 
IStickel et al1ll99lf ). 

Fig. [7] shows the EW distribution used as input for 
our simulation and the resulting one for FSRQs after the 
dilution due to the non-thermal and host galaxy compo- 
nents. The significant shift to lower values is clearly seen. 
This is similar to the comparison between real FSRQs 
and radio quiet QSOs shown in Fig. [3] 

Fig. [8] compares the distributions of f^^, the syn- 
chrotron peak energy, of sources classified as FSRQs and 
BL Lacs in our simu lation with tho s e of b lazars included 
the radio sample of iGiommi et al] (|201ll ). which is the 
sample with the best determination of fp eak values cur- 
rently available. The agreement is clearly quite good 
and reproduces well the fact that BL Lacs tend to have 
^peak values significantly higher than FSRQs. 

We note that, although our input parameters were 
partly based on the WMAP5 sample, it is important to 
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Figure 7. The Mg II EW distribution used as input for our 
simulation (solid line) and the output, diluted distribution for 
FSRQs (dashed line) in the simulated radio survey. 
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Figure 8. Top panel: the fp eak distribution of radio selected 
FSRQs taken from the work of iGiommi et alj l l201lfl (solid 
histogram) compared to that of FSRQs in a simulation of a 
radio survey (dashed histogram). Bottom panel : the i^,,, dis- 
tribution of the BL Lacs in the radio sample of lGiommi et all 
ll201ir ) (solid histogram) compared to that of BL Lacs in a 
simulation of a radio survey (dashed histogram). 



stress that, starting from a single LF and evolution, plus 
a fraction of non-evolving sources, we are able to repro- 
duce quite well the main properties of the two blazar sub- 
classes including those not part of the input (like relative 
fractions and percentage of BL Lacs with redshift). 



5.2 X-ray flux limited survey 

Table [2] summarizes our main results. About 2/3 of our 
blazars are classified as BL Lacs, which is consistent with 
the value of 731^% in the EMSS sample. As in the radio 
case, a small fraction (15%) of the simulated blazars are 
misclassified as radio galaxies. 

The mean redshifts for our simulated FSRQs and BL 
Lacs are in reasonable agreement with the EMSS blazar 
sample values ~ 1 and ~ 0.37. We note that, unlike the 
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1 2 3 

Redshift of featureless BL Lacs 

Figure 9. The redshift distribution of the BL Lacs that show 
a featureless spectrum in our simulation of a radio flux density 
limited survey, and that in a real survey would have no redshift 
determination. 



Table 2. Results from a simulation of an X-ray flux limited 
survey (5 X 10 -13 erg cm -2 s _1 ) 



Source type 




Number of 
sources 


<*} 


(V/Vm) 


FSRQs 




2,836 


1.23 


0.65 


BL Lacs (all) 




5,622 (4,460) 


0.36 


0.51 


BL Lacs (log fp eak 


> 16.5) 


927 (895) 


0.33 


0.45 


BL Lacs (log u s , 

° peak 


> 17) 


185 (177) 


0.34 


0.34 


Radio galaxies 




1,542 


0.04 


0.48 


Total 




10,000 


0.58 


0.55 



WMAP5 sample, the EMSS sample is relatively small (56 
sources) and therefore a detailed comparison is hampered 
by the small number statistics. 

79% of our BL Lacs have a redshift determination, 
in good agreement with the EMSS value of 93+2i%- Al- 
though we assumed that all non-evolving sources do not 
have a standard accretion disk, 30% of the BL Lacs pos- 
sess one and are classified as BL Lacs only because their 
emission lines are swamped by the non-thermal contin- 
uum. The smaller fraction of X-ray selected BL Lacs 
with disks in our simulations, as compared to radio- 
selected ones, is in accordance with the fact that fewer 
EMSS BL Lacs have emission lines clearly detectable 
in their optic al spectra than, f or example, 



1 J y BL 
JckelfeoOll ; 



Lacs (see, e.g. Rector et al. |2000| ; iRector fc Stock c 
IStickel. Fried, fc Kiihr|[l993l ). 

As regards V/Vm, our simulated mean values agree 
reasonably well with the EMSS ones of 0.67 ± 0.08 and 
0.42 ± 0.05 for FSRQs and BL Lacs respectively, deri ved 
using the samples described in lPadovani et all (J2003). 

Fig. [TOl compares the distributions of fp eak of FSRQs 
and BL Lacs in our simulation wi th those of blazars be - 
longing to the soft X-ray sample of lGiommi et alj (|201ll ). 
which includes Planck, Swift and Fermi observed blazars 
and it is therefore probably the sample with the best de- 
termination of z-'peak values currently available. The agree- 
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Figure 10. Top panel: the dist ribution of the X-ray se- 

lected FSRQs in lGiommi et alXll201lh (solid histogram) com- 
pared to that of FSRQs in a simulation of an X-ray survey 
(dashed histogram). Bottom panel: the fp eak distribution of 
the BL Lacs in the X-ray flux limited sample of lGiommi et all 
J201 il l (solid histogram) compared to that of BL Lacs in a 
simulation of an X-ray survey (dashed histogram). 



ment is clearly good reproducing well the fact that BL 
Lacs have much higher ^p eak values than FSRQs. 

5.3 A blazar sequence? 



The existence of a strong anti-correlation be- 
tween bolometric luminosity and fp Ca k, known 
as the "blazar sequence" has been the subject 
of intense d iscuss ions since its first propo sal by 
iFossati et all (Il998t) and iGhisellini et all (Il998l) (e.g. 



Giommi. Menna & Padovani 1 19991 : iPadovani et al 



2003: 
2006: 



Caccianiga & Marcha 2004; Nicppola ct al 



Padovani 



2007 



Ghisellini fc Tavecchid 120081: 



iNieppola et all l2008l : iGiommi et alj 1201 ll) . In this 



tion we use our simulations to comment on the existence 
of such a sequence. 

Figure ITTI shows our radio and X-ray selected simu- 
lated blazars in the log(^p eak ) - log(i/L„(5 GHz)) plane. 
This reproduces the famous plot used by IFossati et al.l 
(1998) to propose the existence of the blazar sequence 
based on the correlation shown in this plane by FSRQs 
and BL Lacs discovered in shallow radio surveys (2 and 
1 Jy samples) and BL Lacs found in the X-ray flux lim- 
ited Einstein slew survey (/ x ^ 10 _12 erg cm -2 s _1 ). In- 
deed, considered together the simulated radio and X-ray 
selected blazars display a broad correlation with radio se- 
lected FSRQs and BL Lacs (red open squares and open 
circles) mostly filling the top left and central part of the 
diagram and X-ray selected BL Lacs (green open circles) 
mostly confined to the lower right corner of the plot. This 
particular positioning of the points (bright FSRQs of the 
LSP type vs. fainter HSP BL Lacs) is not due to any in- 
trinsic correlation between luminosity and f p Cak but is in- 
stead the result of a selection effect resulting from the fact 
that bright radio sources are mostly drawn from the high 
end of the blazar luminosity function, while BL Lacs in 
X-ray flux limited samples are mostly high fp eak sources 
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Figure 11. The FSRQs (open squares) and BL Lacs (circles) of our radio flux density limited (/ r > 0.9 Jy, red symbols) and X-ray 
flux limited (/ x > 5 X 10~ 13 e rg cm" 2 s —1 , green symbols) simulated samples, plotted in the log(y)-log(iAL„) plane, the portion of 
the parameter space used by Fossati et al. (1998) to show the existence of the "blazar sequence" by comparing radio (/ r > 1 — 2 



Jy) and X-ray selected (/ x 



10- 



2 erg cm 2 s 1 ) complete samples. Filled circles represent BL Lacs with very weak lines (EW 



< 2 A) or completely featureless which, in a real survey would most probably not have a measured redshift, and therefore would 
not appear in the plot. Note that almost all of these objects occupy the top right part of the diagram. Blazars in this area are 
high-luminosity HSP sources, that is objects that are "forbidden" in the blazar sequence scenario. 



(intrinsically rare) drawn from the low end of the lumi- 
nosity function where the source density is largest. The 
most imp ortant differe n ce be tween Fig. [TT] and the di- 
agram of iFossatT et all (|l998l ) is in the high-luminosity 
- high fp Cak part, where most of the radio and X-ray 
selected objects with no redshift (red and green filled 
poi nts) are locat e d. Th ese sources could not be plotted 
by iFossati et al.l l|l99ct ) since the luminosity of blazars 
without redshif t cannot be estimat ed. This left the top 
right part of the Fos sati et al.l |T998) diagram empty, thus 
contributing to make the displayed data look like a power 
sequence. 

5.4 Assessing the stability of our results 

To assess the dependence of our results on the adopted 
LF and evolution, we have made the following two checks: 
first, we varied their input values by la adopting $(P) oc 
P" 285 , k = 7.0, and /3 = -1.4 (luminosity peak at 
z ~ 1.9) and <b(P) oc P" 31 , k = 7.7, and /? = -1.7 
(luminosity peak at z ~ 1-75); second, we used as an 
alternative LF the sum of the BL Lac and FSRQ LFs 
based on the beaming model of lUrrv fc Padovanil l|l995l ) 



(converted to Ho = 70 km s _1 Mpc^ 1 ) assuming, for con- 
sistency with the way they were derived, a pure luminos- 
ity evolution of the type P(z) = P(0)exp[T(z)/r], where 
T(z) is the look-back time. A value of r = 0.33, which is 
consistent with the evolution of DXRBS FSRQs and BL 
Lacs combined (based on the samples in iPadovani et al.l 
120071 ), was also assumed. 

We have also checked the dependence of our results 
on the assumed value of the Doppler factor, running 
the simulation with values of (5) ranging from 5 to 20, 
and we have also considered the case of a dependence 
of 5 on radio power as suggest ed by some a uthor s (e.g. 
iHovatta et"aHl2009h . although lEister et all (|2009h . who 
studied a complete sample, have not been able to con- 
firm such a correlation. 

In all cases our main results, that is the prevalence of 
FSRQs and BL Lacs in radio and X-ray selected samples 
respectively, the higher redshifts and evolution of FSRQs 
as compared to BL Lacs, and the no evolution of X-ray 
selected BL Lacs, were confirmed, which shows that they 
are independent of the details of the LF and evolution. 

We also modified our assumption that standard ac- 
cretion disks are associated only with sources having an 
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FR II parent but the simulations we obtained were in- 
consistent with observations. In fact, even if the relatively 
small fraction of 20% of objects with FR I parents are ran- 
domly assigned an accretion disk and allowed to evolve, 
X-ray selected BL Lacs reach too high values of (V/Vm) 
and (z), 0.58 and 0.91 respectively, which are incompat- 
ible with observations. 



6 DISCUSSION 

So far, blazars have been largely discovered in the ra- 
dio, /i-wave, X-ray or 7-ray bands. We have shown that 
the wide diversity of fp ea k observed in blazars results in 
strong differences in the observed flux at radio or X-ray 
frequencies (up to three orders of magnitudes) causing 
severe selection effects in surveys performed in these en- 
ergy bands. This, combined with the fact that blazars 
are always identified on the basis of their optical spec- 
trum, which reflects a mix of non-thermal, thermal, and 
host galaxy emission, has been the origin of a number of 
enduring open issues about the nature of blazars. 



6.1 Evolution 

A long-standing blazar puzzle is the difference in redshift 
distribution and cosmological evolution between FSRQs 
and BL Lacs, selected both in the radio and in the X-ray 
band. BL Lacs are mostly located at low redshifts and ex- 
hibit moderate, or even negative, evolution, while FSRQs 
evolve strongly just like radio quiet QS Os and show a red- 
shift distribution t h at peaks at z > 1 (IStocke et al.lll982l: 
Stickel et ail Il99ll; iRector et al.l 120001: iBeckmann et all 



20031 ; iPadovani et all 120071 : IGiommi et al l 120091 ). Our 

simulations reproduce quite well both of these findings 
(see Fig. [6] and Tables [1] and [2} implying that they are 
due to heavy selection effects. Most of the simulated BL 
Lacs found in radio surveys are luminous objects with 
broad lines that are diluted by non-thermal radiation 
beyond the 5 A EW limit (many of them just below, 
thus allowing a measurement of their redshift), while the 
BL Lacs found in simulated X-ray surveys typically show 
high ^p 0ak values (and therefore are X-ray bright) and 
are drawn from the low-power end of the radio luminos- 
ity function where non-evolving FR Is are preferentially 
found ( Section We note that the (V/Vm) of radio se- 
lected BL Lacs is not too different from that of FSRQs, 
while the {V/V m } of X-ray selected BL Lacs is signifi- 
cantly lower, as found in real surveys (Section 15. ip . 

Another interesting outcome of our simulations is the 
fact that X-ray selected BL Lacs with progressively larger 
values of fp ea k are characterized by lower and lower values 
of (V/Vm) (see Tab. [2J. This is in full agree ment with the 
puzzli ng, and so far unexplained, results of Rector et al.l 
j2000h and of IGiommi. Menna fe Padovanii (|l999T ) who 
reported that the (V/Vm) of BL Lacs is a function of 
their X-ray-to-radio flux ratio. 



of J-'peak energies , which is shifted to values higher tha n 
those of FSRQs ijAbdo et al.ll2010d ; IGiommi et al.ll201ll 1. 
expanding on the well-known fact that high ^p eak objects 
(HSPs) are always BL Lacs. This experimental difference 
is well reproduced in our simulations (see Fig. [S} which 
give (log(^ oak )> = 12.9 for FSRQs and (log(*/ p s cak )) = 
14.1 for BL Lacs for the case of a radio survey. This 
distinction is due to the fact that blazars with higher 
t'pcak values produce more non-thermal optical light than 
low I'peak sources, diluting more easily the broad line 
component, and are therefore classified more frequently 
as BL Lac objects. We note that this has been inter- 
preted in the literature as an intrinsic physical differ- 
ence between LSP (detected mostly in the radio band) 
and HSP (detected mostly in the X-ray and 7-ray band) 
blazars due to the fact that HSPs are observationally 
characterized by a low intrinsic power and external ra- 
diation field, given their very weak or absent emission 
lines. As a consequence, cooling was thought to be less 
dramatic in HSP allowing particles to reach energies high 
enough to produce synchrotron emission well into the X- 
ray band IjGhisellini et all 1998). In our scenario, instead, 
all sources have exactly the same chance of being HSP or 
LSP (that is, the value of 7 pea k is drawn from the distri- 
bution shown in Fig. [4] independently of luminosity) and 
the very different Vp Csii distributions observed in radio 
and X-ray surveys arise from the strong selection effect 
discussed in Section 15.31 and the emission line dilution 
mentioned above. 

Our simulations predict the existence of a significant 
number of BL Lacs with redshift that cannot be mea- 
sured, which occurs when both fp Cak and radio power 
are so large that dilution becomes extreme. For example, 
~ 81% of our simulated sources with ^p eak > 10 15 Hz and 
P T > 10 26 W/Hz have no redshift. This is consistent with 
the fact that most BL Lacs in current 7-ray selected 
samples have no measured redshift, as Fer mi is known 
to preferentia lly select high ^p Cak BL Lacs (|Abdo et al.l 
I2010al . l20lfbl ). This effect is also clearly shown in Fig.lTTI 
where most of the simulated blazars with no measurable 
redshift (filled circles with light colors) occupy the top 
right part of the diagram. Fig. [§] shows the intrinsic red- 
shift distribution of these featureless BL Lacs for the case 
of our simulation of a radio flux density limited survey. 



6.3 What is a BL Lac? 



6.2 



peak 



distribution 



Recent results, based on radio and 7-ray surveys, have 
revealed that BL Lacs, on average, display a distribution 



iBlandford fc Reesl { 19781 ) had originally suggested that 
the absence of broad lines in BL Lacs was due to a very 
bright, Doppler-boosted synchrotron continuum. In the 
years following that paper observations of various BL 
Lacs, mostly selected in the X-ray band, showed that 
in many cases their optical spectrum was not swamped 
by a non-thermal component, as host galaxy features 
were very visible, and it was thought that most BL Lacs 
had intrinsically weak lines. We have shown here that 
these two possibilities are not mutually exclusive and 
indeed are both viable, depending on radio power and, 
therefore, on the band of selection. One important con- 
sequence of our scenario is that objects so far classified 
as BL Lacs on the basis of their observed weak, or unde- 
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tectable, emission lines belong to two physically different 
classes: intrinsically weak lined objects, whose parents 
are LERGs/FR Is (more common in X-ray selected sam- 
ples, since they reach lower radio powers) and heavily di- 
luted broad-lined sources, which are beamed HERGs/FR 
lis (more frequent in radio selected samples). Therefore, 
while the non-thermal engine is probably the same, the 
thermal one is obviously different. This solves at once 
the issue of the FSRQ/BL Lac transition objects and of 
the many differences between BL Lacs selected in the 
radio and X-ray bands, which include line strength, ex- 
tend ed radio emission and morphology, and evolution 
(e.g. iRector fe StockdlioOll . and references therein). This 
hypothesis, which explains also many other open issues 
of blazar research, is directly testable. Our simulations, 
in fact, imply that the majority of sources in high-flux 
density radio-selected samples are identified as BL Lacs 
only because all lines have EW < 5 A in the optical ob- 
serving window. These sources should show a strong (EW 
> 5 A) Ha line if observed in the near or mid-infrared, 
since this is the strongest emission line, and would then 
be considered FSRQs, which reinforces our point on the 
very strong effect that line dilution has on blazar classifi- 
cation. HERG/FR II BL Lacs, then, should be easily rec- 
ognizable through infrared spectroscopy. These sources, 
as expected, are also more dominant at higher radio 
powers: for example, according to our radio simulation 
the fraction of BL Lacs with standard accretion disks is 
only ~ 3% for P r < fO 26 W/Hz but becomes ~ 98% 
above this value. This also means that fainter radio- 
selected samples of BL Lacs should be more and more 
similar to X-ray selected ones, apart from their i/p Cak val- 
ues, since higher ^ ^ are easier to detect in the X-rays 
l|Padovani fc Giommj[l995h . 

The implications of our hypothesis for unified 
schemes is quite straightforward: the parent population 
of BL Lacs need to include both LERGs/FR Is and 
HERGs/FR lis, while that of FSRQs is made up of 
HERGs/FR lis only. This should have only a small effect, 
for ex ample, on the LF fitting done bv lUrrv fc Padovanil 
(|l995h . as HERG/FR II BL Lacs make up the high-power 
end of the radio LF while most of the number density 
comes from LERG/FR I BL Lacs. 



6.4 Blazar classification 

Our new scenario has strong implications on blazar classi- 
fication as well. If the relevant physical distinction for ra- 
dio sources is between LERGs (mostly FR Is) and HERGs 
(FR lis), for most purposes then (LFs, evolution, etc.) 
HERG/FR II BL Lacs should be simply grouped with 
FSRQs. How does one distinguish in practice HERG/FR 
II BL Lacs from LERG/FR I BL Lacs? This is simple in 
the presence of any (even weak) broad lines or for tran- 
sition objects. In other cases (e.g. completely featureless 
spectrum or in presence of absorption features) there is 
no easy way to distinguish between the two subclasses 
although, for example, radio power and/or morphology 
could help. Given the paucity of known LERG/FR Is 
at relative high redshifts, X-ray selected (and also fainter 
radio-selected) BL Lac samples are also useful in selecting 



such sources, which are very relevant also for the study 
of the so-called "AGN feedback" and the role that AGN 
radio emission plays in galaxy evolution through the 
so-called "radio-mode" accretion (|Croton et al.ll2006l ). It 
should also be clear that BL Lacs can be used to study 
the broader issue of the relationship between LERGs and 
HERGs, including their evolution. 

A by-product of our simulations has also been the re- 
alization that some sources classified as radio-galaxies do 
not have their jets oriented at large angles with respect 
to the line of sight, as expected, but are instead mod- 
erately beamed blazars with their non-thermal emission 
swamped by the galaxy (note that none of these objects 
has a standard accretion disk). These sources, which are 
all local (~ 90% at z ^ 0.07) should be recognizable 
by their blazar-like SEDs an d indeed some of them have 



alread y been identified by lDennett- r 
2000 h . lGiommi et al.1 J2002, 2005) and 



2005). 



horoc & Marcha 



Anton fc Browne 



Recently iGhisellini et al.l l|201lf ) have proposed a 
new classification scheme to divide BL Lacs from FS- 
RQs, which is based on the broad line region (BLR) lu- 
minosity in Eddington units and set at a dividing value 
of Lshn/L^dd ~5x 10~ 4 . This turns out to be also the 
value, which separates radiatively efficient (i.e., standard 
accretion disks) from radiatively inefficient (i.e., ADAFs) 
regimes, and therefore coincides with our HERG/FR II - 
LERG/FR I division. Therefore, IGhisellini et~atl <|201ll ) 
are also suggesting that HERG/FR II BL Lacs belong 
with the FSRQs. 



7 CONCLUSIONS 

We have tackled the open issue of the relationship be- 
tween the various blazar subclasses, and of the con- 
sequences of selection effects, through extensive Monte 
Carlo simulations. Our approach is based on robust 
observational input, and on various results obtained 
by many authors over the past 20 years or so, which 
had never before been put together in a comprehen- 
sive way. We kept the number of assumptions and pa- 
rameters to a minimum so as to draw robust conclu- 
sions. Our starting point are two populations of high- 
excitation, high radio-power, evolving and low-excitation, 
low-power, non-evolving radio sources. To these we add: 
an homogeneous synchrotron self-Compton component 
with a distribution of electron Lorentz factors peaked 
at relatively low values, an optical/UV quasar template, 
the EW distributions of the main broad lines observed in 
radio-quiet quasars, a distribution of disk/jet ratios, and 
finally an elliptical host galaxy. 

Our main results can be summarized as follows: 

(I) we explain the main properties of all blazar sub- 
classes, namely FSRQs and BL Lacs selected in the ra- 
dio and X-ray band. Our simulations reproduce well the 
prevalence of FSRQs and BL Lacs in radio and X-ray 
selected samples respectively, the higher redshifts, evolu- 
tion and lower i^p eak of FSRQs as compared to BL Lacs, 
the non evolution of X-ray selected BL Lacs, and the main 
differences between radio and X-ray-selected BL Lacs in 
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a manner that does not dependent significantly on the 
adopted LF and evolutionary details; 

(2) objects classified until now as BL Lacs on the ba- 
sis of their observed weak, or undetectable, emission lines 
actually belong to two physically different classes: most 
BL Lacs in high flux density, radio-selected samples are 
actually beamed radio quasars with their emission lines 
heavily diluted by the non-thermal continuum and for all 
purposes should then be grouped with FSRQs. Infrared 
spectroscopy should reveal in many such sources the pres- 
ence of strong (EW > 5 A) Ha; most BL Lacs selected in 
the X-ray band are instead intrinsically weak-lined, low- 
excitation radio galaxies with a strong non-thermal, jet 
component; 

(3) there are only two main intrinsic blazar types: low- 
ionization (mostly beamed FR Is) and high-ionization 
(beamed FR lis) ones. All other classifications and classes 
proposed so far are not physically relevant and are sim- 
ply due to severe selection effects and different relative 
strengths of the non-thermal, thermal, and galaxy com- 
ponents, which make up blazars' optical emission; 

(4) some sources classified as radio-galaxies on the 
basis of their optical properties are instead moderately 
beamed blazars (that is, sources with their jets forming 
a small angle with respect to the line of sight) with their 
non-thermal emission swamped by the host galaxy. Some 
of these sources have already been recognized from their 
SEDs; 

(5) our simulations show that the purported correla- 
tion between luminosity and fp eak , the so-called "blazar 
sequence" , is likely a selection effect resulting from com- 
paring shallow radio surveys with shallow X-ray surveys. 
We also show that blazars with featureless optical spec- 
tra, and therefore without a redshift determination, are 
mostly high luminosity - high ^p Cak sources, a type of 
blazar that should not exist in the "blazar sequence" sce- 
nario. This is consistent with the lack of redshift in the 
majority of BL Lacs in current 7-ray selected samples, 
as Fermi is known to preferentially select high fp Cak BL 
Lacs; 

(6) the topics addressed in this paper are not relevant 
only to blazars but are also related to the broader issues 
of low-ionization radio galaxies and radio-mode "AGN 
feedback" . 



In this paper we limited our simulations to the radio 
and the X-ray bands where SSC is a fair approximation 
of the observed non-thermal emission. The properties of 
7-ray detected blazars are instead not co nsistent with 
simple SSC models (e.g. lAbdo et~alll2010d l. and almost 
half of the ra dio and X-ray selec ted LSP blazars are 7- 
ray quiet (e.g. lGiommi et al]|201ll ). The present approach 
must therefore be integrated with additional information 
about the properties of the inverse Compton emission 
before it can be used to simulate 7-ray surveys. We are 
planning to extend our simulations to the 7-ray band by 
taking into account the recent results of iGiommi et all 
(|201ll ) who determined the 7-ray properties of blazar 
samples selected in different bands. 
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